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Alternative abstract 

This paper is an invited comment on larXiv:1110.5527l presented at 
Hypercomplex Seminar 2012 and on sixteen earlier published papers 
by Zhidong Zhang and Norman H. March. All these works derive from 
an erroneous solution of the three-dimensional Ising model published in 
2007. A self-contained detailed rigorous proof is presented that the final 
expressions in this work are wrong and that the conjectures on which 
they are based consequently fail. Further errors and shortcomings in the 
follow-up works are also pointed out. 
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Summary 

Thirteen follow-up papers by Zhang and March perpetuate the errors of a 2007 paper by 
Zhang, which was based on an incorrect application of the Jordan-Wigner transformation 
and presents final expressions that contradict rigorously established exact results. The 
presentation given here can be used as a brief mathematical introduction to the Ising 
model for nonexperts. 

1. Introduction 

In a very long paper [Tj published in 2007 Z.-D. Zhang claims to present the exact 
solution of the free energy per site and of the spontaneous magnetization of the 
three-dimensional Ising model in the thermodynamic limit. This claim has been 
shown to be false [2HZ] and we shall show here that very little original work, if any, 
in [1] can be salvaged. 

The principal reason why the outcomes of [1] are wrong is that they contradict 
exactly known series expansion results OH]. Several references were cited in 0E] 
which show that p] violates rigorously established theorems. As these cited theorems 
are formulated for very general lattice models with rather general interactions, re- 
quiring complicated notations and such concepts as Banach spaces and Banach al- 
gebras, it takes some effort to check that every needed detail is there to make the 
proof rigorous. 

Therefore, we present here a simpler self-contained presentation, restricted to the 
three-dimensional Ising model on a simple cubic lattice, which can be used as a short 
introduction for nonexperts interested in this model. 

Definition 1.1. The isotropic Ising model on Z 3 , a periodic n x n x n lattice with 
N = n 3 sites i = {i x , i y , i z ) on a 3-torus, is defined by its configuration space 

(1) T 3 D Z 3 — > {±1} W , i h+ ^ = ±1, for i G Z^, 
and its interaction energy 

(2) n N ■. z 3 -»c, n N = n N ({ai}) = -Jj2 a i°3 - B J2 a *> 

where the sum over is over all nearest-neighbor pairs of sites i and j, J is the 
interaction strength and B is the scaled magnetic field. Sites i and j are nearest- 
neighbor (nn) sites, if and only if 

(3) (i x -3x,iy-jy,iz-jz) = (±1,0,0), (0,±1,0), or (0, 0, ±1) mod n. 
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Remark 1.2. The generalization to the orthorhombic lattice is straightforward, 
replacing n by n, n' , n" and J by J, J', J" for the three lattice directions. We 
consider the isotropic lattice for the sake of simplicity of arguments, as this 
special case suffices to disprove Zhang's claims [T]. 



Definition 1.3. Given a function A = A({a}) of the spin configuration, its 
expectation value is 

(4) {A) N = ±- £ Ae-^, (A)=lim{A) N , 



where the partition function, 

(5) Z N = 



{<Ti=±l} 



is a state sum taken over all 2 N spin configurations, while j3 = (kT)^ 1 with T the 
absolute temperature and k Boltzmann's constant. If (3, J, and B are real, then 
PiW}) = e~@ nN i 'Zn is the Boltzmann-Gibbs canonical probability distribution. 

Definition 1.4. The free energy per site Jn and its infinite system limit / are given 
by 

(6) - 0f N = — log Z N , f = lim /at, 
whereas the spontaneous magnetization is defined by 

(7) / = lim lim (tjj n ) M = lim lim — — 



----- \ v en / /v —— - - „ , is i n vj 

BlON^oo BlON^ooZ N 

{<Ti=±l} 



with iq any of the N lattice sites, as the lattice is chosen periodic. The pair-correlation 
function of spins at sites i and j is (<7j<7j), 



Remark 1.5. As in [T], we shall concentrate on the zero-field (B = 0) thermo- 
dynamic limit (limAr^oo). The order of limits in ([7|) was used implicitly in Yang's 
paper [8] on the spontaneous magnetization of the square-lattice Ising model [9HTT] . 
With the opposite order of limits the result is identically zero. An alternative 
definition is I 2 = lim(<7j<7j) \ B=0 in the limit of infinite separation of sites i and 
3 M- 

In pQ Zhang starts out mimicking the treatment of the two-dimensional Ising mo- 
del by Onsager and Kaufman [9l41 1| . in order to calculate the free energy, 
magnetization and pair correlation of the three-dimensional case. Even though Zhang 
made two early errors in [T] , while transforming to Clifford algebra operators and tre- 
ating boundary terms [5], he claims [6j that these are overcome by two 
conjectures. But these conjectures are based on no serious evidence whatsoever and 
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the resulting expressions for the free energy and magnetization pQ are demonstrably 
incorrect, as they fail the series test [2J[5]. 

First, in section 2, a detailed account will be given of the rigorous results of the 
1960s violated by Zhang's work. Theorems 12.81 and 12 . 91 provide rigorous proof of the 
correctness of the series test. Then, in section 3, further comments will be presented, 
including several on the follow-up work by March and Zhang [3 25 . which contain 
several additional errors and misleading statements. 

2. Some rigorous results of the 1960s revisited 

In recognizing the criticisms to which his work in [1] has been subjected in 
[2j[4l[5l[7] , Zhang (supported more recently by Norman H. March) has argued that 
the usual high-temperature series expansions [26) , renormalization group treatments 
[271128] , and Monte Carlo simulations [291130] , fail to apply in the vicinity of infinite 
temperature owing to singular behavior and Yang-Lee zeros [31][32] present even in 
the thermodynamic limit. Hence, it is argued, such criticisms are not applicable as 
a basis for criticizing the quite different conclusions he has reached. See specifically 
the claims Zhang has made in the second paragraph of [3] , and in the second half of 
page 766 of [6], as well in section 5 of [16] . second half of page 534. The aim of this 
section is to show specifically by a detailed mathematical analysis that there is no 
credibility at all in these claims. 

In fact, five decades ago several theorems were published and supported by fully 
rigorous proofs that underpin the validity of the criticisms of Zhang's work, see, 
e.g., 33,34 for review. Nevertheless, let us here take the reader through a simplified 
treatment especially tailored to apply to the point at issue, namely, the statistical 
mechanics of the Ising model on a cubic lattice with periodic boundary conditions. 

The proof of the thermodynamic (infinite system-size) limit of the free energy 
typically uses the following lemma, see e.g. (2.15) in [33] : 

Lemma 2.1. 

(8) | logTre' 4 — logTre s | < \\A - B\\, for A and B hermitian. 

Proof. The proof follows immediately working out 

logTre- - logTre* = jf' J- logTr e^^dk = | ^gg^ ^, 

where the last integrand is an expectation value. (In this paper we only need to 
consider the commuting case that A and B are diagonal matrices.) ■ 

Theorem 2.2. The free energy per site fx converges uniformly to a limit f as the 
system size becomes infinite for (3 J real and bounded. 



[3] 



Proof. In order to prove this we must estimate — fu | for N, M > No, with 
N sufficiently large. Here we do that only for periodic cubic lattices N = n 3 , 
AI = m 3 and compare with the larger periodic cubic lattice of size NM = (nm) 3 . 
By changing a subset of the interactions we can change the larger lattice into N 
identical cubes of the size M lattice, or the other way around. The proof is then 
provided by counting the changed interactions and by using Lemma [8] In our case, 
the trace in the lemma is just the sum over spin configurations and the norm the 
maximum over all configurations. ■ 

Remark 2.3. Lemma EOl can also be used to show that the free energy / does not 
depend on boundary conditions in the large system limit with different shapes than 
cubes, provided it is taken in the sense of van Hove, see e.g. for details. 



Remark 2.4. The proof of Theorem 12.21 gives a rigorous bound on the difference of 
the free energy per site of a finite system and its large-system limit. It can therefore 
be used to estimate the accuracy of finite-size calculations using e.g. Monte Carlo 
simulations. 



Lemma 2.5. The partition function Zjy ^) is a Laurent polynomial in e@ J , so that 
(3 ffif is singular only for the zeros of this Laurent polynomial and for e@ J — oo. As 
Zn is a sum of positive terms for real (3 J, it cannot have zeros on the real axis. 

We will show that the zero closest to (3 J = (or e' 3,7 = 1) in the complex (3 J plane 
is uniformly bounded away, i.e. Zm ^ for all \(3J\ < Kq and all N for some fixed 
Kq. This means that /at can be expanded in a power series in (3 J that is absolutely 
convergent for \(3J\ < Kq and uniform in N. It is well known that more and more 
coefficients become independent of N as N increases. Together this implies that the 
limiting / also has a power series in (3 J with radius of convergence at least Kq. 

We continue by deriving a lower estimate for Kq. Most proofs of the analyticity of 
free energies and correlation functions use linear correlation identities of Schwinger- 
Dyson type, known under such names as the BBGKY hierarchy, Mayer-Montroll or 
Kirkwood-Salzburg equations. We could use [35] and [36]. But instead, let me give 
an alternative proof using an identity of Suzuki 37,38 , restricted to the isotropic 
Ising model on a simple cubic lattice with periodic boundary conditions and of 
arbitrary size, as this method also can be used to generate the coefficients of the high- 
temperature series. More precisely, using the canonical definition of the expectation 
value of a function A = A({a}) of the spin configuration, we have the correlation 
identity [371155] : 
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Lemma 2.6. (M. Suzuki, 1965 J57$58j) 

(9) (fk) =~i;((iWW* j £>)) . 

\ i=l / AT m fc= i \ V ,=i / V j nn ^ / / N 

where j% , . . . , j m are the labels of m spins and I runs through the labels of the six 
spins that are nearest neighbors of aj k . 

Proof. The proof of |9|) is easy summing over spin aj k in the numerator of the 
expectation value, i.e., 

(10) £ aJ J ^^ l =t^(pj^A E e^-*"*". 

<rj h =±l ^ inn « ' <r jk =±\ 

Averaging over k has been added in ((9]), so that all spins are treated equally, consi- 
stent with the periodic boundary conditions. The lemma is also valid without that. ■ 

Next we use 



Lemma 2.7. 



(11) tghf/3J^cr;j = ai^ai + a 3 ^2a h o- h ai 3 + a 5 ^2a h ai 2 a l3 a h ai 5 , 

^ 1=1 ' (6) (20) (6) 

where the sums are over the 6, 20, or 6 choices of choosing 1, 3, or 5 spins from the 
given a\ , . . . , a&. It is easy to check that the coefficients are 

t(l + 16t 2 + 46t 4 + 16t 6 +i 8 ) -2t 3 



ai (l + i 2 )(l + 6t 2 +t 4 )(l + 14t 2 + t 4 )' ° 3 (l + i 2 )(l + 14< 2 + t 4 )' 
16t 5 

(12) a5= (l+t 2 )(l + 6f 2 +t 4 )(l + 14t 2 + t 4 ) ' tEtgh(W - 

The poles of the a.; are at t — ±i, t = ±(\/2± l)i, and t — ±(\/3 ± 2)i. It can also be 
verified, e.g. expanding the a.; in partial fractions, that the series expansions of the 
ai in terms of the odd powers of t alternate in sign and converge absolutely as long 
as \0J\ < arctg(2 - \/3) = n/12. 

Proof. Clearly, the tanh in (jTTJ) can be expanded as done. Replacing all six 
spins, a i by —ai, shows that no terms with an even number of spins occur. Also, 
permutation symmetry allows only three different coefficients. Multiplying (|ll|l with 
one, three, or five spins oi and then summing over all 2 6 = 64 spin states, is one 
way to derive (|12j) . It is then straightforward to verify the following partial fraction 
expansions, 

1 / Pit P2t \ a/2 / p 3 t p 4 t 

Ql ' 5 ~ 24 V 1 + [ Pl tf + 1 + (p 2 i) 2 J 8 V 1 + (p 3 <)2 + i + (p 4t )S 
, 1 Pht 



(13) a 3 



31+ ( P5 tr ' 

1 / Pit P2t \ 1 P5t 

24Vl + ( Pl <)2 + l + {p 2 tf) 6 1 + (p 5 i) 2 ' 
[5] 



(14) Pi,2 = 2 ± V3, _p 3j4 = V2 ± 1, p 5 = 1, (pip 2 = P3P4 = !)• 
The remaining statements of the lemma follow from these expansions. ■ 



We can now prove the following two theorems for magnetic field B = 0: 

Theorem 2.8. The correlation functions (Y[iLi°~ji)N ond their thermodynamic 
limits (Y[i=i°~ji) are analytic, having series expansions in t or (3 J with radius of 
convergence bounded below by |7?J) and uniformly convergent for all N including 
N = oo . Let d be the largest edge of the minimal parallelepiped containing all sites 
ji > ■ ■ ■ j jm ■ Then the coefficient of t k with k < n— d for the lattice with N — n 3, sites 
equals the corresponding coefficient for larger N, including the one for N = oo. 

Proof. We can assume that m > and even, since for m odd we have (Ili = i o~j t )N 
= as it both is invariant and changes sign under the spin inversion <jj — * — o~i for 
all sites i. 

The system of equations (l9|)~ ([T2l can be viewed as a linear operator on the vector 
space of linear combinations of all correlation functions of the 3-dimensional Ising 
model. It is easy to estimate the norm of this operator. Using the alternating sign 
property of the cij's, it is easy to verify that a\, a%, and a§ can all be written as 
t times a series in i 2 , which three series consist of positive terms only when t is 
imaginary. This means that each |aj| is maximal for given \t\ when t is imaginary 
and within the radius of convergence, i.e. p2 in (|14j) . 

From the 32m terms in the right-hand side (RHS) of © after applying it 
follows then that we only need to study 

for purely imaginary t to find the desired upper bound r for the norm. Setting t = \x 
with < x < 2 — \/3 to stay within the first pole of (fT5|) , we next define 

™ 2-^(3 -x 2 )(l-3x 2 ) 

(Id) r = t, ^rr^ 7"r^ 7T; tor X = \t\. 

y J (1 - x 2 ){\ - Ux 2 + x 4 ) 11 

We then have that the RHS of ([9]) is bounded by rM, where M = max \(a ■ ■ ■ a) \ with 
the maximum taken over all 32n pair correlations in the RHS. (Obviously, M < 1 if 
(3 > and real, but we shall not use this.) We can easily show that r < 1 for 

|i| < (V3 - \/2)(\/2 - 1) = 0.131652497 • • • , or 
(17) \PJ\ < arctg[(\/3 - V2)(V2 - 1)] = 0.130899693 

To prove analyticity of (JlHi a ii)N as a function of /3 at (3 = 0, we apply ^ to it. 
Then we apply ([9]) to each of the 32m new correlations, and we keep repeating this 
process ad infinitum. Since a 2 = 1, we will from some point on regularly encounter 
the correlation with m = 0, i.e. zero a factors, for which (1) = 1, so that the iteration 
process ends there. Each other correlation (with m > 0) vanishes with at least one 
power of t, as can be seen comparing e.g. (|9]) and (fT2|) . We conclude that the iteration 
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process generates the high-temperature power series in t to higher and higher orders, 
for arbitrary given size N of the system. 

To get the partial sum of the series to a given order, we only need to keep the 
contributions for which the iteration process has ended and expand all occurring etj 
as series in t. The sum of the absolute values of the terms is bounded by ^2 rJ < oo 
when (|17|) holds. However, the original correlation function is meromorphic with a 
finite number of poles away from the real t axis for any finite N. Thus for sufficiently 
high order of series expansion in t, the remainder term is arbitrarily small. The only 
possible conclusion is that we have proved convergence of the series expansion of 
(nl=i a ji)N in powers of t, uniform in N with a finite radius of convergence in the 
complex t and (3 planes bounded below by (fT7]) . 

To prove the final statement of the theorem for finite N, we notice that the above 
iteration process generates new correlations with the range of the positions j of the 
spins extended by one in a given direction. As long as we do not go around a cycle 
(periodic boundary condition) of the 3-torus, we do not notice any TV-dependence. 
It takes at least n — d iteration steps to notice the finite size of the lattice. 

Combining the convergence uniform in N with the fact that more and more 
coefficients converge with increasing N, we conclude that (Q™i a ii)N converges 
to a unique limit as N — > oo for \t\ < 2 — V3 5 with the properties stated in the 
theorem. ■ 



Theorem 2.9. The reduced free energy /3/jv for arbitrary N and its thermodyna- 
mic limit (3f are analytic in (33 for sufficiently high temperatures. They have series 
expansions in t or (3J with radius of convergence bounded below by [lTjj and 
uniformly convergent for all N including N — oo. The first n — 1 coefficients of 
these series for N — n 3 equal their limiting values for N — oo. 

Proof. To prove analyticity of j3f in terms of f3 at f3 = it suffices to study the 
internal energy per site or the nearest-neighbor pair correlation function, as 

(18) un = —(Hn) n = ~^7p~ = _ 3J(cto,o,oo'i,o,o)jv! 

as follows from ([5]) and ©■ Here 00,0,0 an d 01,0,0 can be any other pair of neighboring 
spins. The proof then follows from Theorem 12.81 and integrating the series for un, 
using Z N \p =0 = 2 N , implying lim^o /3/jv = - log 2. ■ 



Remark 2.10. Adding a small magnetic field B and generalizing the steps in the 
above, we can conclude that all correlation functions are finite for small enough \/3\ 
and \j3H\, so that there are no Yang-Lee zeros j31[l32] near the H = axis for small 
(3 and H. The proof can also be generalized to the case that the interactions are 
anisotropic, i.e. J, J', J" as in [T]. Then Zn is a Laurent polynomial in each of e@ J , 
e 13 - 7 , and e^ J , etc. 
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Remark 2.11. Similar results can be derived for the low-temperature series, for 
example after applying the Kramers- Wannier duality transform to the high-temper- 
ature regime of the dual system with spins in the centers of the original cubes and 
with four-spin interactions around all cube faces perpendicular to the edges of the 
original lattice [39] . 

Remark 2.12. It is straightforward to calculate the first few high-temperature 
series coefficients of the free energy by the method described in this section, with or 
without using the averaging in (|9|) . They agree with the long series reported in [26] 
and earlier works cited there. Zhang's free energy formula claimed for all finite 
temperatures p] does not agree, as already the coefficients of k 2 = t 2 in (A12) and 
(A13) of p] differ. Zhang's excuse that there are two expansions, one for finite (3 and 
one for infinitesimal f3, violates general theorems, that apply to more general models 
than the Ising model [IK]. Here this excuse is invalidated in detail by Theorem l2.9l 

Remark 2.13. Zhang's spontaneous magnetization series is obviously wrong. In 
three dimensions one should have I — 1 = 0(x 6 ), with x = e~ l3J in the low- 
temperature limit, x — > (J > 0), as each spin has six nearest neighbors [5J, rather 
than eight, which would result in the four-dimensional I — 1 = 0(x 8 ) presented by 
Zhang in (103) of pQ. 

Remark 2.14. The finite radius of convergence of the series expansions about /3 = 
is also hinted at by the fact that the zeros of Zn for T = oo occur for B = ±ioo, 
[3B = ±m/2 [7J. For fixed temperature T or j3 = 1/kT and [3. J > real the zeros 
of Zn lie on the unit circle in the complex cr 2 ^ B plane 32J, all located at — 1 
at infinite temperature [7J and spreading out with decreasing temperature until 
the zeros "pinch" +1 on both sides of the unit circle at and below the critical 
temperature, in agreement with the theory of Yang and Lee [3U|32]. Zhang's claim 
that this pinching at +1 also occurs at = [HElElIlB] is disproved by Theorem 

m 

Corollary 2.15. As pointed out already in '2,4,51, all final results of 0/ are proven 
wrong, as they do not agree within a finite radius of convergence with the well- 
known series expansion coefficients. This also means that the conjectures of are 
falsified. 
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3. Further remarks and objections 

3.1. Two series expansions for the same object 

In appendix A of pQ Zhang claims to reproduce the first 22 terms of the high- 
temperature series for the free energy. But this is no more than reverse engineering, 
fitting the known coefficients [26 to an integral transform (A.l) or (74) in [1\ giving 
the first few coefficients of the weight functions as given in (A. 2). There is no more 
information than the series results provided by others, so that this does not constitute 
a new result, as explained in [U[S]. 

As this construction this way is based on a conjectured integral transform of 
weight functions that can only be reconstructed from a few known series coefficients, 
it cannot be considered an exact solution. Knowing this, Zhang conjectures ad hoc 
above (A. 3) on page 5400 another choice for the weight functions, namely w x = 
1, w y = w z = 0, leading to another high-temperature series for non- infinitesimal 
temperatures, in violation of the rigorous result on the uniqueness of the series 
expansion presented in section 2. This is not sound mathematics [2]. 



3.2. Citations by other authors 

The outcomes of pQ have been criticized in [29 , 3D] , as they disagree with recent 
high-precision Monte Carlo calculations presented there. Both the position of the 
critical point and the values of the critical exponents differ from the ones in [1], 
while the results of [29,30] agree with those of many others obtained by a variety of 
methods [28] . 

One paper on a decorated three-dimensional Ising model [40] , mapping this model 
exactly to the Ising model on a cubic lattice, used Zhang's free energy [1] as an 
approximate result in the analysis. The experimental paper |41j states that their 
result for the critical exponent A = 2.0 ± 0.5 is consistent with pQ. However, the 
reported error bar is so large that this means nothing. Moreover, paper [42] on the 
Heisenberg model only briefly cites [1] as an Ising reference. 

The authors of [43] [44] learned from pQ the quaternion setup of the transfer 
matrix of the three-dimensional Ising model, which was well-known earlier, see e.g. 
[45H47] HI In Zhang's work pQ this is treated before the first error occurs with the 
Jordan- Wigner transformation to Clifford algebra operators. His P's and Q's do not 
anticommute [5j[6]. 



1 Maddox only published his final formula for the free energy [45]. The details were 
discussed in a special session, where also his error (the same as Zhang's first error [5]) was 
discovered. 
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3.3. Advertising wrong critical exponents 



Klein and March took the exact Ising critical exponents for dimensions d = 1,2,4 
together with the proposal of [T] for d — 3 and made an ad hoc fit [IS] for all real 
1 < d < 4. However, they failed to compare with the results from e-expansion [2"Tll2"8] . 
where e = A — d. This is a serious shortcoming, as the [48] formulae disagree with the 
e-expansion exponents for small e and fail the one foremost explicit test available. It 
may also be noted that the extrapolated Ising exponents for d = 3 from e-expansion 
agree with those extracted from series expansions and Monte Carlo calculations [55] , 
while differing from those presented in pQ . 

March and Zhang have followed this paper [48] up with thirteen publications, 
thus perpetuating the errors of the original work pQ. Some of these works compare 
Zhang's critical exponents with those from experiments on CrBr3 and Ni [131114] . 
Nickel is known to have Heisenberg exchange interactions and its critical exponent 
(3 is about the accepted value for the three-dimensional Heisenberg model, which is 
also about Zhang's value wrongly claimed for Ising. 

Comparing experiments with models needs a discussion of the interactions in 
the experimental compounds, whether Ising or Heisenberg, isotropic or anisotropic, 
short-range or long-range, etc. No such analysis was presented. The same objection 
can be brought up about section 2 of [20] . 

In [15] critical exponents for the two- and three-dimensional g-state Potts model 
are discussed. Those for d = 2 are by now well established, but the values presented 
for d — 3 cannot all be correct, as for the Ising case q — 2 the exponents of [1] have 
been used. 

In [18, 19J a new formula for critical exponent <5 is given, improving the one in [48] . 
The same objection still applies, as again no comparison with e-expansion is made. 

It is implied by the theory of Yang and Lee [SUES], that the best experimen- 
tal results on Ising exponents are to be expected from measurements on liquid-gas 
transitions in simple substances. March and Zhang have admitted that the expo- 
nents of pQ fail this test, see section 3 of [20]. Their suggestion that the experiment 
needs to be redone carries no credibility, as the critical exponents measured in a 
number of similar experiments are indeed typical for Ising, see section 3.2.2 of [28] . 

3.4. Singularity of free energy at T — oo 

Several statements in section 5 of [16] repeat and expand on statements in [T1G2E] 
contradicting rigorous theorems discussed in section 2 above. For example, while it 
is correct that the free energy / diverges at T = oo, this does not correspond to a 
physical singularity, as the combination (3f is to be used. Indeed, e~@f relates to the 
normalization of the Gibbs ensemble probability distribution and j3f is the principal 
object of Theorem 12.91 Multiplying fjf with kT results in / having a convergent 
Laurent expansion with a leading pole at T = oo that has no physical significance. 
Another point is discussed in Remark l2.14l in section 2. 
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3.5. False argument for a. — 



Paper |T7| addresses tricritical behavior. The authors claim that the logarithmic 
divergence of the specific heat, a = (log) at tricritical points in three dimensions, 
supports the similar value reported in [1] for the Ising critical behavior. However, 
this reasoning is flawed lacking any theoretical basis and contradicts the accepted 
value a = 0.110 ± 0.001, see eq. (3.2) and tables 3-7 of [28]. 

3.6. The e — d — 2 expansion 

In papers ;21,[52] on Anderson localization the authors say that e = d — 2 is not 
a small parameter for d = 3, just like e = 4 — d of the e-expansion is not. This 
ignores that the best e-expansion extrapolation results agree remarkably well with 
those from series, Monte Carlo, and experiment [2S]. This argument to support pQ 
is again not valid. 

3.7. Higher dimensions 

The combinatorial sums defining the 3-dimcnsional Ising model involve commuting 
spin variables and an interaction energy that is a function of these spins. There is 
no reason to introduce time and quantum mechanics in this classical system, as is 
done in |23j . On the other hand, introducing the transfer matrix changes one space 
coordinate to (discrete) imaginary time. After "Wick rotation" to real time the 3- 
dimensional Ising model relates to a (2+l)-dimensional quantum system. The fourth 
dimension introduced in [T] has only been used to obtain wrong results violating 
rigorous results. 

3.8. Fractal dimensions based on wrong results 

In [24 Zhang and March write down some proposals for fractal dimensions. However, 
the values given for dimension 3 are based on incorrect results of pQ. 

3.9. Unfounded Virasoro algebra 

The most recent paper [25] uses the weight factors of p] to introduce a Virasoro 
algebra in 3+1 dimensions. This is ad hoc and the notations in equation (5) and 
seven lines below (6) there are not mathematically sound. To take the real part 
of the absolute value of a phase factor instead of just writing 1 makes no sense. 
Also, the Virasoro algebra relates to an infinite dimensional symmetry, which is only 
consistent with conformal symmetry in two dimensions, see e.g. [49] . Therefore, [25] 
has fundamental errors. 
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Open letter to Dr. Zhang in response to arXiv:1209.3247 



After posting arXiv:1209.3247 Dr. Zhang sent me an email ending with "Any further 
criticisms or comments from you are welcome." The following is a I^TgX version of 
the body of my September 17 email reply: 

First of all, the comment arXiv:1209.0731 contains very little that is in my publi- 
shed earlier Comment and Rejoinder. There is some material from the unpublished 
additions to the arXiv version of the Rejoinder, but that part is much improved with 
several new details added. There are also some pages discussing papers published 
later. Therefore, the statement "hard to find anything new" is wrong. Also, several 
statements are not even addressed in the Response and cannot be covered with the 
"unnecessary to repeat all of Zhang's responses" . 

Conjecture 1 is not backed up by any quantitative evidence in the original 117 
page work. The validity can at this moment only be judged by the resulting free 
energy. In my section 3.1, I make clear that the result becomes an integral transform 
from the weight functions to the free energy, but that there is no clear argument to 
get the weight functions. You obtained two series based on two choices, one fitting 
series (A2) to the high-temperature series to as many terms as known, the other 
choosing (1,0,0) leading to a different series. The first way you have no result, as 
you have no more than the known series terms. 

The second way can be disproved by series, providing the first few terms of 
the well-known high-temperature series are rigorously established. The older 
proofs are correct but not easy to read. Therefore, I gave a much simpler proof 
with mathematical precision. 

My proof does not depend on the papers by Lebowitz and Penrose and by 
Gallavotti et a!., contrary to what you seem to suggest. Also, the statements about 
these papers are taken out of context: 

The statement Re/3 > on page 102 in Commun. Math. Phys. 11, 99 (1968) is 
needed when the gas model has no hard core. Section II opens with the statement 
that analyticity at (3 — can be shown for a hard core potential. The Ising model 
is equivalent with a lattice gas version with at most one particle per lattice site 
(empty-occupied becomes spin +/ — 1), a special case of a hard core on the lattice. 
Thus this objection of yours does not apply. 

The statement about (3 > in [10] on page 3 of your response arXiv:1209.3247 
does not appear explicitly in [10]. It appears near the bottom of the left column of 
page 494 of Phys. Lett. A 25, 493 (1967). This is a very short letter without much 
detail. To prove a finite radius of convergence you need to prove an inequality with 
a positive j3. Then (3 — will be included within the radius of convergence. Your 
statement means that the Gallavotti et al. letter was poorly written at that point, 
[probably due to printer error as Phys. Lett, did not allow authors to correct proofs], 
not that it is wrong. Also, as said before, I used nothing of these two papers, so that 
these objections are misplaced. 
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The next objection that (3f does not exist at T — oo is also invalid. The combi- 
nation (3f — — In 2 there, as / has a simple pole at (3 — 0. One finds that / has a 
Laurent expansion with pole term (— ln2)//3 followed by a power series in (3 with a 
finite radius of convergence. Statistically, at (3 = all states have equal probability 
and there is no phase transition, as not only /?/, but also all correlation functions 
are analytic at (3 — 0, in spite of the fact that interactions are turned on. 

When you expand A = Z X I N in (A12) and (A13) of the original 117 page 
paper, you expand exp(— f3f), which is equivalent to expanding (3f. This makes your 
objection to expanding (3f unreasonable. The finite radius of convergence proof that 
I gave proves that (A13) is not correct. 

You bring up that 1/Z has a zero at (3 = in the infinite system. But this is 
again misleading. Yang-Lee theory is only about the zeroes of partition function Z: 
When zeros pinch the real temperature axis in the large system limit, then there is 
a phase transition. There is no theorem for 1/Z. 

This pinching cannot occur, as the proof given by me can be extended to the 
double expansion of /?/ in /? and f3B. The proof for the more general cases is in the 
old literature. From this joint analyticity at /3 = and (3B = 0, it follows that zeros 
are a finite (nonzero) distance away, contradicting the pinching that you claim. 

Remark 2.4 in my comment and similar statements now make it possible to 
test your free energy with Monte Carlo methods, as one can now estimate both 
the systematic error due to finite size and the statistical error due to Monte Carlo 
methods. 

Also, the latest various experimental and theoretical estimates for a are signifi- 
cantly different from 0, see the review of Pelissetto and Vicari, which is [28] in my 
comment. Please, check the accuracies reported of large numbers of theoretical and 
experimental works that are discussed there. 

Next, the use of dimensionless Ki = (3Ji and h = (3H can be done in more than 
one way. The partition function and correlation functions (and (3f) only depend 
on these combinations. That some authors set (3 = 1, does not mean a loss of the 
high-temperature case. If you have the result in the i*Q and h, you also can choose 
a new /3, say (3' , and write in the results Ki = (3 1 Ji and h — (3'H. There is no loss 
of the high-T limit (3' — 0. Again this is an objection that is invalid and it does not 
apply to my latest comment, as I nowhere used (3 = 1, nor did I use results from 
authors that did. 

Finally, the last paper arXiv:1110.5527 is based on an incorrect solution of the 
3D Ising model. There are problems I noted: With <p a phase, \e l ^\ = 1, thus formula 
(4) contains phases that drop out. Also, having three independent Virasoro alge- 
bras means that you have the 3+1 dimensional space rewritten as a 6-dimensional 
(2+2+2) space, the "product" of three independent 2-dimensional spaces. Things 
do not add up. 

This is my answer to your "Any further criticisms or comments from you are 
welcome." 
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